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Abstract - We consider a single spin-down impurity atom interacting via an attractive, short- 
range potential with a spin-up Fermi sea in two dimensions (2D). Similarly to 3D, we show how 
the impurity can form a metastable state (the "repulsive polaron" ) with energy greater than that 
of the non-interacting impurity. Moreover, we find that the repulsive polaron can acquire a finite 
momentum for sufficiently weak attractive interactions. Even though the energy of the repulsive 
polaron can become sizeable, we argue that saturated ferromagnetism is unfavorable in 2D because 
of the polaron's finite lifetime and small quasiparticle weight. 



^ [ Introduction. — The Fermi gas with effectively re- 
I pulsive, short-range interactions provides a model system 
■ O m which to investigate magnetic instabilities of the Fermi 
^ liquid state. Of particular interest is the putative Stoner 
^transition to itinerant ferromagnetism, which is driven 
I— purely by the repulsive interactions. Such a scenario 
can potentially be realized and investigated with ultra- 
j^eold, dilute gases of atoms [1]. However, recent cold-atom 
experiments in three dimensions (3D) have cast doubt 
T-H ©n whether the Stoner transition can ever be realized in 
\l a strongly repulsive Fermi gas with short-range interac- 
^"^tions P]. The central issue is that strongly repulsive in- 
f eractions in 3D can only be truly short-ranged if the un- 
derlying potential is attractive. This necessarily demands 
^^that any magnetic phase be metastable at best, since the 
• . true ground state will involve attractively-interacting pairs 
_^6f t and i fermions. In 3D, the ferromagnetic state is ap- 
parently never stable [2^4] , thus a key question is whether 
this is also true in lower dimensions. 
. P.' 

Lowering the dimensionality will strongly modify both 
the two-body scattering properties and the many-body 
properties. In ID, it is in fact permissible to have purely 
repulsive contact interactions, but here the Lieb-Mattis 
theorem precludes any ferromagnetic ground state [5]. 
The situation is less clear in 2D, however, where one 
still requires an underlying attractive potential to generate 
strong, short-range repulsive interactions, but in contrast 
to 3D there is always a two-body bound state and the 
two-body scattering amplitude is energy dependent, even 



at low energies [B]. 

Fortunately, 2D Fermi gases have recently been realized 
in cold-atom experiments [7l413|. making it an ideal time 
in which to investigate ferromagnetism in 2D. Theories 
based on a mean-field approach currently support the ex- 
istence of itinerant ferromagnetism |14j . However, given 
the importance of quantum fluctuations |14| . we ideally 
want to explore limits of the problem where the theory 
is more controlled. To this end, we consider the limit 
of extreme spin imbalance, where we have one spin-down 
"impurity" immersed in a Fermi sea of spin-up particles. 
Here, theoretical approaches based on simple variational 
wave functions |151I16| have proved to be extremely reli- 
able in 3D when compared with experiments |171I18[ and 
quantum Monte Carlo simulations [1111201 ■ Thus far, the 
ground state of the 2D impurity problem has been theoret- 
ically investigated piH^ , and it has even been suggested 
that recent 2D experiments have already observed single- 
impurity physics |24j . 

In this Letter, we show how a J,-impurity interacting at- 
tractively with a t Fermi sea in 2D can form a metastable 
state — the so-called "repulsive polaron" — with energy 
greater than that of the non-interacting impurity. Our fo- 
cus is on equal masses (m-f = wj,), given its relevance to 
current cold-atom experiments and ferromagnetism, but 
we also investigate how the repulsive polaron depends on 
the mass ratio mi/mj,. As the attraction is decreased, we 
find that the energy of the repulsive polaron increases 
and eventually surpasses the Fermi energy ej? of the spin- 
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up Fermi sea. Moreover, we find that the repulsive polaron 
acquires a finite momentum for sufficiently weak attrac- 
tion. However, its decay to the ground-state "attractive" 
polaron is also enhanced when — sp, which suggests 
that saturated fcrromagnctism is unlikely to exist in 2D, 
as we argue below. 

Methods. — In the following, we consider a two- 
component (t, 2D Fermi system with short-range at- 
tractive interactions, described by the Hamiltonian 
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where — /2m„ (with h = 1), Q, is the system area 
and g is the bare coupling constant describing the inter- 
species contact interaction. In addition we note that, for 
fermions, the interaction is interspecies only since the ex- 
clusion principle forbids intraspecies s-wave scattering. In 
2D, the bare coupling strength can be related to the two- 
body binding energy eg, which is always present for at- 
tractive interactions in 2D, via [B] 
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The cut-off momentum A can be sent to infinity at the 
end of the calculation. 

To proceed, we adopt the variational polaron wavefunc- 
tion introduced in Ref. 1151. 



IP) 



at\l,\FS) 



'^kq''cp+q_k;Ck-fCq-f 

kq 



\FS). (3) 



The first term on the r.h.s. describes a bare impurity and 
the second term takes into account one particlc-hole pair 
excitation of the Fermi sea, \FS). The polaron ground- 
state energy is determined by minimising {P\H\P) while 
keeping {P\P) ~ 1. The variational approach can also 
be generalized to study metastable excited states such as 
the repulsive polaron by instead minimising the action 
/ dt{P\ (ihdt-H) \P) — see, e.g., Ref. [55] for further de- 
tails. Equivalently, one can use the diagrammatic method 
in which the interspecies interaction is treated in the lad- 
der approximation [25J . This yields the expression for the 
self-energy: 
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ejj and E measures the energy 
change due to the addition of an impurity to the Fermi 
sea. Here, and in what follows, the hole and particle mo- 
menta are limited to |q| < kp and |k| > kp, respectively, 
with kp the majority Fermi momentum. In principle, one 
can obtain a more accurate description by allowing more 
particle-hole pair excitations; however, it has been shown 



that one particle-hole pair excitation yields a good esti- 
mate in 3D, as we will discuss below. 

From the self-energy (U), one can obtain the energy of 
the polaron by locating the quasiparticle pole. 
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Solutions of this equation are in general complex, the pres- 
ence of an imaginary part indicating a finite lifetime of the 
quasiparticle. If the decay rate F is much smaller than 
9^e{e}, the quasiparticle energy is given by 



E = 
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We quantify the momentum dependence of the energy E 
with the effective mass m* = l/[d'^E/dp'^]p=o. The quasi- 
particle weight is described by the residue 
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The approximation becomes exact in the limit T/E 0. 

Equation ([5]) yields two solutions for the energy E. The 
negative energy solution E- corresponds to the attractive 
polaron studied in Refs. [2T1[22]. For a weakly interact- 
ing system, the attractive polaron is the ground state of 
an impurity immersed in the Fermi sea. However, as the 
strength of the attractive interaction is increased, a bound 
diatomic molecule dressed by particle-hole excitations be- 
comes energetically favorable [21]. Additionally, a solu- 
tion of Eq. ^ with positive energy E^ exists as found in 
Ref. [23] . The theory of this repulsive polaron was recently 
elucidated in 3D [?7H2^ . This metastable quasiparticle is 
dressed by a cloud of particle-hole pair excitations whose 
effect is, in contrast to the attractive reduced ma- 

jority particle density around the impurity [28] . However, 
the repulsive polaron is not the true ground state (indeed, 
it is not generally an eigenstate of the system) and thus 
it will eventually decay. Energy and momentum conserva- 
tion restrict the possible decay modes: The simplest are 
into either an attractive polaron and a particle-hole pair 
or into a molecule, two holes and one particle [TOll^l^ . 
the first being the main source of decay in the region of 
interest as we shall argue below. 

The decay rate of the repulsive polaron sets the 
timescale at which the physics related to this excited state 
may be observed. In general, if the quasiparticle self- 
energy is known exactly, the rate may be calculated as 
F = — Z+CJm{S(p, However, when the self-energy 
is approximate, a reliable estimate requires the lower ly- 
ing states into which the quasiparticle can decay to be 
incorporated correctly. The self-energy obtained from the 
wavefunction ^ is unable to accurately capture the decay 
to an attractive polaron (described by the same wavefunc- 
tion) and a particle-hole pair. Instead, we calculate the 
rate of decay to the attractive polaron following the idea 
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Fig. 1: (Color online) Energy of the repulsive polaron as a func- 
tion of eb/sf for different experimentally relevant mass ratios 
(thick): From top to bottom, m^/m_i = 40/6, 1, 6/40. Where 
these curves split, the dotted lines depict the zero-momentum 
energy while the lower branch traces the minimum of the po- 
laron dispersion. We also display the asymptotic result pO[l for 
equal masses (double-dash-dotted). At the bottom, the thin 
(black) lines are the energies of an infinitely massive impurity 
as obtained from the exact solution (solid) and the ladder ap- 
proximation (dotted). As eb/ef — >■ 0, both of these approach 
ef- 



of Ref. [2H]: Close to the quasiparticle poles, the repulsive 
(-I-) and attractive (— ) polarons are described by Green's 
functions 
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The decay to the attractive polaron is dominated by pro- 
cesses which correctly match the energy difference IS.E = 
E+ — E-. Thus, we estimate the decay rate by replacing 
the bare impurity propagator in the ladder summation by 
the pole expansion of the attractive polaron. In this man- 
ner we arrive at the expression 



F = -Z+3mY^ 
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The energy E^^.^ is obtained from -Ekqip by replacing 
with the effective mass of the attractive polaron, and the 
couplingconstant g is adjusted to obtain ultraviolet con- 
vergence!. 

Repulsive polaron. — The properties of the repulsive 
polaron are parameterized by the dimensionless quantity 
Eb/sf, which gives a measure of the interaction strength. 



^We expect the impurity to act as a free particle for momentum 
k ^ kF- Instead of using a complicated dispersion for the attrac- 
tive polaron (cfTectively letting depend on momentum), we note 
that the large k behavior is cut off by a corresponding term in the 
coupling g and adjust the coefficient in front of this term to cancel 
the ultraviolet divergence. 



Fig. 2: (Color online) Inverse effective mass of the repulsive 
polaron for mass ratios m^/m_|_ =40/6 (dot-dashed), f (solid), 
6/40 (dashed). Where m\ < 0, a polaron with finite momen- 
tum will be energetically favorable. Ref. [24] considered the 
equal mass case in the regime eb/ef ^ 3 and their results 
agree with ours. Inset: (rn^ = m_|_) Momentum P* which min- 
imizes the energy of the polaron. 



For strong attraction, eb/sf ^ 1, the repulsive polaron 
branch is only weakly perturbed by the two-body bound 
state which sits far below the continuum, and the repul- 
sive polaron approaches the bare impurity state. Con- 
sequently, when Eb/ef oo, the energy — !■ 0, the 
residue — > 1 and the effective mass — >■ toj, , as 
depicted in Figs.[l][3J respectively. Indeed, in this limit we 
find 

E+/ef ^ {m^ + m;)/[TO; log{EB/SF)] (10) 

to leading order, which is consistent with previous per- 
turbative calculations for equal masses [30]. The two- 
body decay rate F is also logarithmically suppressed for 

Et » ef: From Eq. (®, F ^ ""l^rT'X ,^^^^7%-' 

" ^1—1" m-^lm-f+ml) [log{eB/eF)r 

which approaches zero faster than E^ as expected (note 
that in this limit Z_ oc Ef/sb)- In this limit, we find that 
the rate TpM at which the repulsive polaron can form a 
molecule is even smaller, like in 3D [28]. Conservation laws 
require the creation of an additional particle-hole pair, and 
Fermi antisymmetry leads to a suppression by a further 
Ef/eb by the same argument as in 3D [21] . Thus the rate 
of the resulting three-body process is expected to approach 
zero at least as fast as Vpm/£f ~ {ef/ebY- Eventually, 



when the attractive polaron is no longer the ground state 
(e.g. above Eb/Ep ~ 10 for equal masses [H]), the pole 
expansion ((HJ) breaks down for the attractive polaron and 
the repulsive polaron will predominantly decay into the 
molecular ground state. 

For decreasing eb/ef: the energy Ej^ initially steadily 
increases until it eventually exceeds Ef (see Fig. [T]). The 
interaction at which this occurs depends sensitively on 
the mass ratio; For equal masses, the critical interac- 
tion {eb/£f)c — 2.7, while for m^lm:^ oo we have 
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{£b/£f)c oo, and for m^/mi = we have {£b/£f)c = 
0. Note that our resuh in the latter case is consistent with 
the exact solution E+ = J^" ^ cot^^ [i log ^] and 
in fact we see in Fig. [T] that is quite close to the exact 
energy, apart from where the exact energy is non-analytic 
near Sb/^f ~ 0. Indeed, we expect our variational ap- 
proximation ([3]) to be least accurate for the case of an 
infinitely massive impurity, where higher order particle- 
hole excitations of the Fermi sea are expected to be more 
important [16) . 

For weak attraction, Sb/sf 0, we find that the ef- 
fective mass of the repulsive polaron diverges (Fig. [5]), 
signalling that the lowest-energy repulsive polaron has a 
finite momentum P* beyond that point. This was also 
concluded from the shape of the spectral function calcu- 
lated in Ref. [23]. From the inset of Fig. [5J we see that 
P* smoothly evolves from to kp as Eb/^f ^ 0- This 
behavior is tied to the fact that, in this limit, the residue 
decreases towards zero and the wave function evolves 
into a bare impurity at momentum p — kp plus one 'f- 
particle excited from zero momentum to kp. This limit 
mirrors the repulsive-polaron solution for infinite impu- 
rity mass, where one removes a f-particle from the bound 
state formed from the zero-momentum state and places it 
at the Fermi surface. In terms of the wave function (jSj, 
the final state at Eb/ef ~ corresponds to a'^^ = and 

'^kq' ~ '5q.o<^k,-fci7pi which is orthogonal to the attractive 
polaron (with Z_ = 1). Thus, in this limit, the energy 
is minimised when p ~ kp, giving E-i- = Ep. Contrast 
this with the p = state, where i?+ = Ep{l -\- m^/m^). 
Both these states are eigenstates of the non-interacting 
Hamiltonian and thus the decay rate of the repulsive po- 
laron approaches zero as eb/ep (see Fig. [3]) . Note 
that when m*^ < 0, the p = state may excite a particle 
out of the Fermi sea to gain momentum P* . By evalu- 
ating the available phase space close to the divergence of 
TO^, the rate at which this process occurs is found to go 
as [E+{\p\ = 0) - E+{\p\ = P*)]5/4, which implies that 
the p = state is metastable with respect to the finite 
momentum repulsive polaron in this limit. 

Similar finite-momentum excitations are encountered in 
the corresponding impurity problem in ID. The Bethe 
ansatz solution for attractive interactions and equal 
masses exhibits an excited "broken-pair" state with neg- 
ative effective mass [33]. Here, the excitation's lowest en- 
ergy is always at momentum kp and equals Ep. By con- 
trast, the energy at p = exceeds Ep, with a maximum of 
2£p when the interactions are switched off, like in 2D. The 
Bethe ansatz solutions are, of course, exact eigenstates of 
the ID system and thus the decay rates for these excita- 
tions are zero. However, in 2D, the repulsive polaron is, 
at best, metastable outside of the limits Eb/ep and 
£b/£f oo, and indeed we see in Fig. [3] that the decay 
rates are sizeable in the regime 1 < Eb/ef < 10 for a 
range of mass ratios. 




Fig. 3: (Color online) Decay rate F of the repulsive polaron 
for mass ratios m-^/mi = 40/6 (dot-dashed), 1 (solid), 6/40 
(dashed). Where these lines split, the upper and lower curves 
depict r for the finite- and zero-momentum polaron, respec- 
tively. Inset: (m^ = mj,) Quasiparticle residues. The solid 
line follows the lowest lying repulsive polaron and the dotted 
line the zero-momentum state. Ref. |M] considered the zero- 
momentum states and our results are consistent. 



Itinerant ferromagnetism. — To examine how this 
impacts the existence of itinerant ferromagnetism in 2D, 
we start by assuming the presence of spin-polarized do- 
mains and then assessing the stability of these. In order to 
have stable, well-defined t and ], domains, we require there 
to be an energy cost for transporting particles across the 
'\-\. interface. Thus, wc require a repulsive polaron formed 
from a J, (t) impurity to have energy greater than the 
Fermi energy of the 4- (t) domain. Assuming mechanical 
equilibrium, where the pressures of each domain are equal, 
we then obtain the stability condition Ej^ > E^yJln^Jml 
{E+ > Ep^/m^^/■mf). For equal masses, this simply gives 
E+ > Ep, which corresponds to Eb/ep < 2.7 as described 
earlier. At first sight, this might suggest that ferromag- 
netism can exist in this regime; however we see in Fig. [3] 
that the decay rate F is large near Eb/ep = 2.7 and 
amounts to a significant fraction of the Fermi energy. The 
large F corresponds to a large uncertainty in the energy 
of the repulsive polaron, thus allowing particles to tunnel 
across the interface and destabilize the domains even when 
£'+ > Ep. Once across, the particles quickly decay into at- 
tractive polarons, e.g. for the typical densities used in 2D 
experiments |10| . a repulsive polaron with T ^ O.Iep near 
sb/sf ~ 1-3 only has a lifetime of order 1ms. Of course, 
the decay rate is eventually suppressed once Eb/ep 0, 
leading one to suspect that ferromagnetism becomes possi- 
ble in this limit. However, is also heavily suppressed as 
Z_ — !> 1, which means that the particles will mostly tunnel 
directly into the attractive polaron state (the probability 
of tunneling into each polaron state is proportional to Z±). 
Thus, our results indicate that saturated ferromagnetism 
is highly unstable in 2D. 
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Another route to determining the existence of itinerant 
ferromagnetism is to start with a uniform 50/50 mixture of 
spins and then compare the formation rate of spin domains 
with the rate of decay into t-4- pairs/dimers. In 3D, pre- 
vious theoretical studies [21131 have shown that the decay 
into pairs always dominates, thus ruling out the sponta- 
neous formation of ferromagnetic domains. Our stability 
argument, however, places a stronger restriction on ferro- 
magnetism in 2D, namely that saturated ferromagnetism 
cannot exist even if fully-polarized spin domains were to 
be artificially engineered. 

A similar scenario holds for unequal masses. In Figs.[T][3] 
we show the properties of a ^''K [^Li] impurity immersed 
in a ^Li [^°K] Fermi sea. In this case, we find the do- 
main stability condition to be satisfied for es/eF.Li < 5.7. 
However, again we find that the decay rate is prohibitively 
large and phase separation is not favored. Our results are 
derived under the assumption that the system is correctly 
described by the Hamiltonian ([1]). However, Feshbach res- 
onances in this heteronuclear system are narrow |34| which 
changes the two-particle scattering properties. Studying 
this problem in detail is beyond the scope of this letter, but 
the enhanced closed-channel character close to a narrow 
resonance decreases the two-body decay rate of the repul- 
sive polaron [35] . While in general the energy will also be 
lower, there may be an intermediate regime of parameter 
space which favors phase separation. 

Our arguments may be extended to realistic experimen- 
tal conditions where the atoms arc subjected to a weak 
trapping potential in the 2D plane. The stability of the 
spin polarized domains relies on mechanical equilibrium 
at the t-i interface. This constraint is local and thus the 
above considerations are also valid in the trapped system 
provided the trap changes smoothly enough that the local 
density approximation holds. Additionally, the atoms are 
confined to quasi-2D by a strong transverse field character- 
ized by the frequency lOz ■ For the equal- mass case studied 
in Ref. jlOj with experimental parameters = 27r x 80 
kHz and e^? = 27r x 9 kHz, we have es ^ Ef ^ ijJz in the 
regime of interest for the question of possible ferromag- 
netism. Thus, we conclude that only the lowest transverse 
mode is occupied and the 2D approximation considered in 
this paper is valid. 

Concluding remarks. — We now turn to the ques- 
tion of the accuracy of the variational wavefunction ([5]). 
In 3D, the energy and effective mass of the ground-state 
quasiparticle have been measured in Refs. [T7] and [15] . 
respectively, and good agreement with the variational cal- 
culation [TSJ[25] was obtained, even in the limit of strong 
interactions. The approach has been further validated 
by agreement with Monte Carlo simulations [20], while 
Ref. [16] has argued that the accuracy of the wavefunc- 
tion arises from a nearly perfect destructive interference 
of the contributions of states dressed by more than one 
particle-hole pair. This latter argument does not depend 
on dimension and indeed it can be shown that the ground 



state energy from the variational approach agrees well 
with the exact Bethe Ansatz solution in ID [55]. While the 
variational wavefunction allows an accurate calculation of 
ground state properties of the system, it is, perhaps, less 
intuitive that this should be the case for a metastable po- 
laron. Nevertheless, recent experimental and theoretical 
studies [3S] of a ""^K-^Li mixture in 3D have demonstrated 
that not only is the energy and residue of the repulsive po- 
laron correctly captured in the equivalent diagrammatic 
approach, the decay rate is also well approximated by a 
technique similar to the one which led to Eq. ©• 

To conclude, we have calculated the quasiparticle prop- 
erties of the repulsive polaron. We find that in the limit 
of weak attractive interactions the lowest lying repulsive 
polaron has finite momentum. These properties may be 
observed by radiofrequency spectroscopy |24j . However, 
the fast decay and small residue of these quasiparticles in 
the regime of strong repulsion indicates that spin domains 
will be unstable and precludes itinerant ferromagnetism in 
the same manner as in the 3D experiment [2] . Our results 
are also applicable to a bosonic impurity immersed in a 
Fermi sea and thus have implications for phase separation 
in 2D Bose-Fermi mixtures. 
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